The media effect in Axelrod's model explained 
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We revisit the problem of introducing an external global field - the mass media - in Axelrod's 
model of social dynamics, where in addition to their nearest neighbors, the agents can interact with 
a virtual neighbor whose cultural features are fixed from the outset. The finding that this apparently 
homogenizing field actually increases the cultural diversity has been considered a puzzle since the 
phenomenon was first reported more than a decade ago. Here we offer a simple explanation for it, 
which is based on the pedestrian observation that Axelrod's model exhibits more cultural diversity, 
i.e., more distinct cultural domains, when the agents are allowed to interact solely with the media 
field than when they can interact with their neighbors as well. In this perspective, it is the local 
homogenizing interactions that work towards making the absorbing configurations less fragmented 
as compared with the extreme situation in which the agents interact with the media only. 
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PACS numbers: 87.23.Ge, 89.75.Fb, 05.50.+q 



I. INTRODUCTION 

Although computational social science - a branch of 
computer science that aims at predicting the behavior 
of groups as well as exploring ways of influencing that 
behavior - is currently in vogue in the scientific media 
(see, e.g., [J), it has a long record in the computational 
physics community 0-13] ■ In fact, physicists have offered 
a variety of models of social influence and opinion for- 
mation (see, e.g., but they seem to have adopted 
as a paradigm a model proposed by a political scientist 
- Axelrod's model @ - which exhibits the right balance 
between simplicity and realism, and so succeeds in boil- 
ing down a collective social phenomenon to its functional 
essence fTol|. In addition, Axelrod's model has served as 
inspiration to the proposal of heuristics that use social in- 
teractions to solve combinatorial optimization problems 

Axelrod's model was introduced to explore the mech- 
anisms behind the persistence of cultural and opinion 
differences in a society, given that the interactions be- 
tween people tend, on the average, to make them more 
alike in their beliefs and attitudes. Hinging on a few 
reasonable assumptions, Axelrod's model offers a non- 
trivial answer to that issue. In that model, an agent is 
represented by a string of cultural features, where each 
feature can adopt a certain number of distinct traits. 
The interaction between any two agents takes place with 
probability proportional to their cultural similarity, i.e., 
proportional to the number of traits they have in com- 
mon. The result of such interaction is the increase of 
the similarity between the two agents, as one of them 
modifies a previously distinct trait to match that of its 
partner. In spite of this homogenizing assumption, the 
model exhibits global polarization, i.e., a stable multicul- 
tural regime Q. Subsequent analysis of Axelrod's model 
by the statistical physics community has revealed a rich 
dynamic behavior with a nonequilibrium phase transi- 
tion separating the global polarization regime from the 
homogeneous regime, where a single culture dominates 



the entire population [T^ - [l6| . 

Another counterintuitive result derived from the anal- 
ysis of Axelrod's model is that the introduction of an 
external global media (i.e., it is the same for all agents) 
aiming at influencing the agents opinions actually favors 
polarization [17| . This result conflicts with the view that 
mass media, such as newspapers and television, are in- 
struments to control people's opinions and so homogenize 
society [l^. The effect of the mass media in Axelrod's 
model has been extensively investigated (see, e.g., [T9I - 
\2^). but no first-principles explanation for it was offered. 
In fact, much of that research has focused on the search 
for a threshold on the intensity of the media influence 
such that above that threshold, the population is cultur- 
ally polarized and below it, the population is homoge- 
neous. More recently, it was argued that this threshold 
is an artifact of finite lattices, and that even a vanishingly 
small media influence is sufficient to produce cultural po- 
larization in a region of the parameter space where the 
homogeneous regime is dominant in the absence of the 
media 23] . 

Here we offer a simple explanation for the media effect 
in Axelrod's model, which is based on viewing the prob- 
lem from a new perspective: rather than focusing on the 
counterintuitive domain-fragmenting effect of increasing 
the media strength, we centre at the trite homogenizing 
effect of the interactions between agents, which become 
more important as the media strength is weakened. We 
conclude then that it is the local homogenizing interac- 
tions that work towards making the absorbing configu- 
rations less heterogeneous as compared with the extreme 
situation in which the agents can interact with the media 
only. 



II. MODEL 

In the original Axelrod's model Q, each agent is char- 
acterized by a set of F cultural features which can take 
on q distinct values. The agents are fixed in the sites of 
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a square lattice of linear size L with free boundary con- 
ditions (i.e., agents in the corners of the lattice interact 
with two neighbors, agents in the sides with three, and 
agents in the bulk with four nearest neighbors). These 
were the boundary conditions used in the seminal paper 
by Axelrod [sj , which we adopt here because they greatly 
facilitate the implementation of Hoshen and Kopelman 
algorithm for counting the number of clusters in a lat- 
tice [131 • The initial configuration is completely random 
with the features of each agent given by random inte- 
gers drawn uniformly between 1 and q. At each time we 
pick an agent at random (this is the target agent) as well 
as one of its neighbors. These two agents interact with 
probability equal to their cultural similarity, defined as 
the fraction of common cultural features. An interac- 
tion consists of selecting at random one of the distinct 
features, and making the selected feature of the target 
agent equal to its neighbor's corresponding trait. This 
procedure is repeated until the system is frozen into an 
absorbing configuration. 

The introduction of a global media in the standard 
model follows the ingenious suggestion of adding a vir- 
tual agent which interacts with all agents in the lattice 
and whose cultural traits reflect the media message [l7l |. 
In the original version, each cultural feature of the vir- 
tual agent has the trait which is the most common in 
the population - the consensus opinion. To speed up the 
simulations, here we choose to keep the media message 
fixed from the outset, so it really models some alien in- 
fluence impinging on the population. Explicitly, we gen- 
erate the culture vector of the virtual agent at random 
and keep it fixed during the dynamics 19, 20]. The inter- 
action of the media (virtual agent) with the real agents 
is governed by the control parameter p G [0,1], which 
may be interpreted as a measure of the strength of the 
media influence. As in the original Axelrod's model, we 
begin by choosing a target agent at random, but now it 
can interact with the media with probability p or with 
its neighbors with probability 1 — p. Since we have de- 
fined the media as a virtual agent, the interaction follows 
exactly the same rules as before. The original model is 
recovered by setting p — Q. 

The simulation of the dynamics of Axelrod's model 
can be greatly accelerated by introducing a list of ac- 
tive agents [la]. An active agent is an agent that has 
at least one feature in common and at least one feature 
distinct with at least one of its neighbors, included the 
virtual agent (media). Clearly, since only active agents 
can change their cultural features, it is more efhcient to 
select the target agent randomly from the list of active 
agents rather than from the entire lattice. Note that the 
randomly selected neighbor of the target agent may not 
necessarily be an active agent itself. In the case that the 
cultural features of the target agent are modified by the 
interaction with its neighbor, we need to re-examine the 
active/inactive status of the target agent as well as of 
all its neighbors so as to update the list of active agents. 
The dynamics is frozen when the list of active agents is 



empty. 

Because the cultural features of the virtual agent are 
fixed a priori, there are three distinct types of absorb- 
ing states of the dynamics. These absorbing states de- 
termine the stationary (or frozen) regimes of the model. 
First, there is the homogeneous regime where all agents 
are identical to the media. Second, there are [q - ly dis- 
tinct homogeneous regimes in which all F cultural traits 
of the agents are distinct from those of the media. How- 
ever, we will show below that this type of absorbing con- 
figuration is absent for large lattices. Third, there is a 
heterogeneous regime characterized by the co-existence 
of domains of agents whose cultural traits are identical 
to the media with domains of agents whose traits are all 
distinct from the media. 

As pointed out before, the puzzling phenomenon re- 
ported first by Shibanai et al. [l7j and then confirmed 
by many other authors (l9l - [23j ] is that for a set of con- 
trol parameters, say F = 5 and q = 5, the homogeneous 
regime in the absence of media is replaced by a hetero- 
geneous regime whenever the media is turned on (i.e., 
for p > 0). This goes against the intuitive idea that the 
effect of the media is to make the culture more homoge- 
neous. Studies based on lattices of small sizes (typically 
L = 40) have suggested the existence of a threshold value 
Pc on the intensity of the media influence such that for 
p > Pc the stationary regime is heterogeneous, whereas 
it is homogeneous for p < Pc [l9l - [2^ . (We recall that the 
parameters q and F are chosen such that the population 
is culturally homogeneous for p = 0.) This finding was 
questioned in a recent study which suggested that pc is 
an artifact of finite lattices and hence that even a vanish- 
ingly small media infiucncc is sufficient to destabilize the 
culturally homogeneous regime [23l |. However, the prob- 
lem with that study was that the virtual agent (media) 
did not take part in the determination of the active sta- 
tus of the agents. As we will argue below, correction of 
this point does not change the claim that pc ~> when 
L — > oo. 



III. SIMULATION RESULTS 

We consider first the ratio 5 G [0, 1] between the num- 
ber of clusters or cultural domains S and the lattice 
area (see fig. [T|). A cluster or domain is simply a 
bounded region of uniform culture. Two or more cultural 
domains characterized by the same culture are counted 
separately in the calculation of S. For the parameter set 
F = q = 5, analysis of small lattices indicates a threshold 
value Pc ~ 0.03 [l^HO], but the results exhibited in fig. 
[T] indicate unequivocally that the heterogeneous regime 
displaces the homogeneous one even for very small values 
of the media strength as the lattice size increases. This 
figure is virtually indistinguishable from that obtained in 
the case the media does not take part in the determina- 
tion of the active status of the agents [1^. In this, as 
well as in the next figures of this paper, the error bars 
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FIG. 1: Logarithmic plot of the ratio g between the number 
of domains and the lattice area as function of the linear size 
L of the lattice for (top to bottom) p = 1, 0.05, 0.04, 0.03, 0.02 
and 0.01. The solid straight line is which corresponds 

to the value of g in the uniform regime. The parameters are 
F = 5 and 5 = 5. 

are smaller or at most equal to the symbol sizes. We note 
that there are many more cultural domains in the case 
the local interactions are turned off (i.e., p = 1) than in 
any other situation. 

A more lively illustration of the effect of increasing the 
lattice size is offered in fig. [2l which shows the fraction 
of runs S,h which ended up in a homogeneous absorbing 
configuration. Typically we carried out from 10^ to 10* 
independent runs of the stochastic dynamics for a given 
value of the media strength p. This figure highlights the 
somewhat odd predominance of the homogeneous absorb- 
ing configurations for lattices of intermediate size, which 
has led to the suggestion of a nonzero threshold value 
Pc- For p — 1 and L > 5 we find that the absorbing 
configurations are always heterogeneous. In addition, we 
note that, except for small lattice sizes, the homogeneous 
absorbing configurations are identical to the media. In 
fact, we have estimated the fraction of runs which 
ended up in the media configuration and plotted in fig. 
12] the fraction of homogeneous absorbing configurations 
distinct from the media, i.e., the quantity 1 — £,m/^h- As 
shown in the figure, this fraction vanishes exponentially 
fast with the increase of the lattice size and so, for 
large lattices, there is essentially a single homogeneous 
absorbing configuration - the one where all agents as- 
sume the media cultural traits. 

The results exhibited in figs. [1] and [2] merely reassert 
the claim that the media promotes cultural diversity, 
since the number of cultural domains increases with the 
strength of the media, measured by the parameter p. 
However, this is a rather odd view of the situation. In 
fact, in all runs and for all parameter settings considered, 
we found that the largest domain is composed of agents 
identical to the media. In addition, the size of this do- 



FIG. 2: Fraction of runs trapped into homogeneous absorbing 
configuration as function of the linear size L of the lattice 
for (top to bottom) p = 0.01, 0.02, 0.03, 0.04, 0.05 and 1. The 
parameters are F = 5 and q = 5. 
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FIG. 3: Fraction of the homogeneous absorbing configurations 
which are distinct from the media as function of the lattice 
size for (top to bottom) p = 0.01, 0.02 and 0.03. The 
solid lines are the fittings exp (— fcpL'^) . The parameters 
are F = 5 and q = 5. 



main is larger than the sum of the sizes of the other 
domains, as shown in fig. SI which exhibits the average 
fraction 77 of the lattice sites (agents) that do not belong 
to the largest domain. Clearly, = 1 — (Smax) /L"^ where 
{Smax) is the average size of the largest domain (media) 
which, together with g, has often been used as an order 
parameter for Axelrod's model [l^l- We note that 
although some of the small domains are also character- 
ized by the media traits, they count as distinct domains 
as they are disconnected from the largest cluster - a sit- 
uation known as cultural diaspora [25|- In agreement 
with the conclusion drawn from the previous figures, the 
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FIG. 4: Fraction of the lattice sites that are not part of the 
largest domain (media) as function of the linear size L of 
the lattice for (top to bottom) p = 1, 0.05, 0.04, 0.03, 0.02 and 
0.01. The horizontal line is (1 — l/g)^.The parameters are 
F = 5 and q = 5. 
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FIG. 5: Fraction of the lattice sites that are not part the 
largest domain (media) in the thermodynamic limit, rjoo, as 
function of the media strength p. The solid line is the fitting 
»7oo = 1820p*'^. The parameters are F = 5 and q — 5. 



sum of the sizes of the domains different from the media 
occupy a finite (i.e., macroscopic) portion of the lattice 
in the thermodynamic limit. In this limit, the average 
size of the domains, excluded the largest one, can be ob- 
tained by calculating the ratio rj/g and yields about 3 
sites, regardless of the value of p, so they are microscopic 
domains. This implies that there is only one macroscopic 
domain in the system - the media - which in the limit 
p — > tends continuously to the size of the entire sys- 
tem. This point is shown explicitly in fig. [S] where we 
plot the asymptotic values of the fraction of lattice sites 
that are not part of the largest domain 7700 against the 
media strength p. 

Figure 0] hints that there is something trivial about 
the nature of the cultural domains distinct from the me- 
dia. The are two points which are worth emphasizing 
here. First, when the agents are allowed to interact solely 
with the fixed external media, the absorbing configura- 
tions are heterogeneous (see fig. since agents whose 
F cultural traits are all distinct from the media at the 
onset will be unaffected by the dynamic rules. For in- 
stance, the horizontal line in fig. U] gives the fraction of 
sites distinct from the media regardless of the domain 
organization in the case the local interactions are turned 
off {p — 1), i.e., (1 — . Those sites are frozen from 

the very beginning of the simulation. Second, there are 
many more agents different from the media for p = 1 
than for p < 1, i.e., the local homogenizing interactions 
contribute to increase the number of agents with the me- 
dia traits. Summing up, due to the non-zero similarity 
restriction for the occurrence of interactions, the media 
in Axelrod's model does not promote uniformity; on the 
contrary, there are more cultural diversity when the me- 
dia is the only partner the agents can interact with than 



when local interactions are allowed. 



IV. CONCLUSION 

There is nothing extraordinary about the effect of the 
media in Axelrod's model, provided we interpret the re- 
sults summarized in fig. 2] from top to bottom: when 
assembling the random initial configuration a fraction 
(1 — of the agents have their F cultural traits dis- 

tinct from the media and so they form the cultural do- 
mains in the case the local interactions are neglected 
(p — !■ 1). As the media strength p decreases, the local 
interactions have the chance to flip some of those traits, 
resulting in the increase of the size of the media cluster. 
Finally, when p ^ the media domain takes over the 
entire lattice, leading to a unique homogeneous absorb- 
ing configuration. For p = the homogeneous absorb- 
ing configuration is not unique - there are different 
possibilities - but, as in the ferromagnetic Ising model, 
a vanishingly small external field suffices to break this 
symmetry and direct the relaxation to a unique absorb- 
ing configuration. 

In terms of the order parameters g (the ratio be- 
tween the number of domains and the lattice area ) and 
(Smax) /L'^ (the ratio between the size of the largest clus- 
ter and the lattice area) the scenario is as follows. Except 
for p — g is nonzero and so the number of domains 
grows linearly with the lattice area - this is the media 
induced diversity. However, for all values of p there is a 
single macroscopic cluster - the media cluster - the size 
of which increases as the media strength p decreases and 
ultimately encompasses the entire lattice for p = (ex- 
plaining this counterintuitive finding is the goal of this 
paper). All other domains are of finite size (typically 3 
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sites in the thermodynamic limit) but they are so nu- 
merous that the sum of their sizes grows with (see 
fig. E]). Therefore, within the perspective that there is a 
single macroscopic domain we may say that the system 
is ordered in the presence of the media regardless of the 
value of its strength p. 

Although we have restricted our analysis to a single 
setting of the control parameters, namely, F = q = 5, 
due mainly to the computational demand of the simula- 
tions of large lattices, our conclusion is valid in general 
since it provides an intuitive first-principles explanation 
for the effect of the media in Axelrod's model, which is 
not affected by the tuning of the control parameters of 
the model. 

The somewhat trivial nature of the cultural domains of 
Axelrod's model in presence of a external media is remi- 
niscent of the or gani zation of the domains of the extended 
majority model |26l - l28l |. which is completely determined 
by the initial random assignment of cultural traits. As 
in that model, relaxation to an absorbing configuration 
is very fast since the random clusters formed at the be- 
ginning are practically unaltered by the update rules of 
the model. 

In summary, we show here that the effect of exter- 
nal global media in Axelrod's model can be described 



much more prosaically by pointing out that the local ho- 
mogenizing interactions destabilize the domains of agents 
(usually, isolated sites) which have zero similarity with 
the media and hence are not affected by it. As a re- 
sult, decreasing the media strength and, consequently, 
increasing the role of the local interactions will favor less 
fragmented absorbing configurations. Reversely, increas- 
ing the media strength will favor more heterogeneous ab- 
sorbing configurations. It seems to us that the cause of 
the ruckus about the media effect in Axelrod's model was 
the failure to realize that the absorbing configurations 
can already be very fragmented due to the interactions 
with the media only. In fact, the mere inclusion of the 
data for p = 1 in figs. [T] and S] turns the media effect from 
a surprising phenomenon to a trivial one. 
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